Estimation of the control parameter from symbolic sequences: Unimodal maps with 

variable critical point 
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The work described in this paper can be interpreted as an application of the order patterns of 
symbolic dynamics when dealing with unimodal maps. Specifically, it is shown how Gray codes can 
be used to estimate the probability distribution functions (PDFs) of the order patterns of parametric 
unimodal maps. Furthermore, these PDFs depend on the value of the parameter, what eventually 
provides a handle to estimate the parameter value from symbolic sequences (in form of Gray codes), 
even when the critical point depends on the parameter. 



In this paper, the order patterns of unimodal 
maps are studied. It is shown how to construct 
order patterns of unimodal maps from their sym- 
bolic dynamics with respect to the partition of 
the state space introduced by the critical point. 
Finally, it is shown that for a subclass of para- 
metric unimodal maps, the study of those order 
patterns allows to estimate the parameter of the 
map that has generated the symbolic sequence. 



I. INTRODUCTION 

Sarkovskii's theorem shows that order and dynamics 
are intertwined in one-dimensional intervals. It is there- 
fore not surprising that the study of the ordinal structure 
of deterministic time series gives valuable information on 
the underlying dynamical system. This work focuses on 
the reconstruction of the so-called order patterns of cer- 
tain unimodal maps, from "coarse-grained" orbits in form 
of 0-1 sequences: if the corresponding iterate lies to the 
left of the critical point, and 1 otherwise. Such binary 
sequences will be called Gray codes. The relationship be- 
tween the Gray codes of parametric unimodal maps and 
the value of the parameter that controls a particular dy- 
namic, was shown in [1, 2, 3]. Other important tool for 
the understanding of one-dimensional dynamical systems 
is the study of their order patterns [4]. Indeed, order pat- 
terns allow to distinguish chaos from white noise, and can 
provide useful information on the parameter or param- 
eters controlling the dynamic of chaotic systems. The 
main goal of this paper is to estimate the control param- 
eter of unimodal maps by means of their order patterns 
alone, even when the exact values of their orbits are not 
accessible but only the corresponding Gray codes. 

The rest of the paper is organized as follows. First of 
all, the general framework is set in Sect. II. In Sect. Ill, 
the concept of order pattern is introduced, and its depen- 
dence on the control parameter is analyzed for the logistic 



and the skew tent maps. Sect. IV summarizes the theory 
on Gray codes. How the order patterns of unimodal maps 
are obtained using Gray codes is explained in Sect. V; its 
application to control parameter estimation is explained 
in Sect. VI. The results presented in this paper are re- 
capitulated in Sect. VII, where some final comments are 
also included. 



II. SCENARIO 

The work described in this paper focuses on a class of 
unimodal maps, hereafter denoted as T . A map f : I —* 
I, where / = [a, b] C R, a < b, belongs to the class T if 
it satisfies the following conditions. 



1. / is continuous. 

2. f(a) = f(b) = a. 

3. / reaches its maximum value /„ 
terval [a m ,b m ] C I, a m <b m . 



< b in the subin- 



4- /(/max) < x c , where x c is the middle point of the 
interval \a m ,b m ], i.e., x r ~ 



5. f(x c ) > x c . 

6. / is strictly increasing function on [a, a m ] and 
strictly decreasing on [b m ,b\. 

The class T includes maps defined in a parametric way, 
say, f\(x) = ip(X,x), where x € I = [a, b], A G J C ffi is 
called the control parameter, and ip is a self-map of / x J. 
Two different situations are considered in this paper: 

1. The control parameter determines the maximum 
value of the map. In this case, the parametric func- 
tion fx is given by 



f x (x) = \F(x), 



(1) 
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where F € T and F(x c ) = F max . The subclass of 
maps f\ £ T complying with this description will 
be denoted by T\ . 
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2. The control parameter is the value of the critical A G J C K (as it occurs for /a G J^i, JF 2 ), and the set P„- 
point, i.e., x c = A. This leads to a new subclass of is defined as 
maps Ti- 

Prr = {x G I : x is of type 7r} , (3) 




where 7r G Sl, then depends on f\ and, consequently, 
on A. According to the ergodic theorem [7, p. 34], if f\ is 
crgodic with respect to the invariant measure fi, then the 
orbit of x G / visits the set P^ with relative frequency 
(i (Pjt), for almost all x with respect to \x. As a result, it 
is possible to study the dependence of P v on A by count- 
ing and normalizing the occurrences of tt in sliding win- 
dows of width L along Of x (x), x being a 'typical' initial 
condition. In the following two subsections this is done 
experimentally with the logistic map (as representative 
of T\) and with the skew tent map (as representative of 
T-z). Since we are primarily interested in the relation be- 
tween the probabilities (J,(P W ) (or relative frequencies) of 
order patterns tt G Sl and the control parameter A of the 
map considered, we will refer to it as the A-distribution 
function (in short: A-DF) of tt, since they are related to 



x the probability distribution functions (wc fix tt instead of 

fixing A). 

FIG. 1: fx( x ) for k = 0, 1, 2, 3 and the corresponding order 
patterns of length 4 for the logistic map when A = 4. 



III. ORDER PATTERNS 

Given a closed interval / C R and a map /:/—>/ 
, the orbit of (the initial condition) x G / is defined as 
the set O f (x) = {f n (x) : n G N }, where N o = {0}UN = 
{0,1,...}, f°(x) = x and f n (x) = /(/" _1 (a;)). Orbits 
are used to define order L-patterns (or order patterns 
of length L), which arc permutations of the elements 
{0, 1, L — 1}, L > 2. We write tt = [no, ni, . . . , 7T£_i] 
for the permutation i— » t:q, L — 1 i— ► 7rj,_i. 

Definition 1 (Order pattern). The point x £ I is said 
to define (or realize) the order L-pattern tt = tt(x) = 

[tT , 7Ti, . . . , 7TL_l] if 

r«(x)<n(x)<...<r^( x ). (2) 

Alternatively, x is said to be of type tt. The set of all 
possible order patterns of length L is denoted by Sl- 

For further reference, it is convenient to assign an inte- 
ger number to each order pattern. This can be made, for 
instance, by means of the Trotter- Johnson algorithm [ ] . 
The order patterns of length 4 along with their "ordering 
numbers" , are shown in Table I. 

As emphasized in [6], there always exist order L- 
patterns with sufficiently large L that are not realized 
in any orbit of / G T . These order patterns are called 
forbidden patterns, whereas the rest of order patterns are 
called allowed patterns. In general, if f\ is a family of 
self-maps of the closed interval / Cl parameterized by 




Order patterns zu 3.7 



FIG. 2: Relative frequency of the order patterns realized by 
the logistic map when L — 4 and A G [3.7, 4]. 



A. Order patterns for the logistic map 

The logistic map, defined as 

f x (x) = Xx(l - x), (4) 

for x G [0,1] and A G [1,4], belongs to T\. The logistic 
map with A = 4 was studied in [6, 8] from the ordinal 
point of view. In Fig. 1 the allowed order 4-patterns 
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# 


Order pattern 


# 


Order pattern 


# 


Order pattern 


# 


Order pattern 





[0, 1, 2, 3] 


1 


[0, 1, 3, 2] 


2 


[0, 3, 1, 2] 


3 


[3, 0, 1, 2] 


4 


[3, 0, 2, 1] 


5 


[0, 3, 2, 1] 


6 


[0, 2, 3, 1] 


7 


[0, 2, 1, 3] 


8 


[2, 0, 1, 3] 


9 


[2, 0, 3, 1] 


10 


[2, 3, 0, 1] 


11 


[3, 2, 0, 1] 


12 


[3, 2, 1, 0] 


13 


[2, 3, 1, 0] 


14 


[2, 1, 3, 0] 


15 


[2, 1, 0, 3] 


16 


[1, 2, 0, 3] 


17 


[1, 2, 3, 0] 


18 


[1, 3, 2, 0] 


19 


[3, 1, 2, 0] 


20 


[3, 1, 0, 2] 


21 


[1, 3, 0, 2] 


22 


[1, 0, 3, 2] 


23 


[1, 0, 2, 3] 



TABLE I: Order patterns of length four. 



for the logistic map with A = 4 arc shown. For this 
value of the control parameter there exist twelve allowed 
order patterns. However, the main goal of this paper is to 
analyze the relationship between the control parameter of 
maps in T\ or T2, and their order patterns, what calls for 
the distributions of allowed patterns for different values of 
A. Figure 2 depicts the relative frequencies of each order 
4-pattcrn for A £ [3.7,4], the patterns being labeled as 
in Table I. To be more specific, for every A, a sufficiently 
long orbit was generated, the occurrences of the different 
order patterns were counted using a sliding window of 
width 4, and finally the counts obtained were normalized 
by the number of windows. These results are estimates of 
the probabilities for the corresponding order patterns to 
occur. Let us point out that, since the physical invariant 
measure of the logistic map is only known for A = 4, 
numerical estimation of those probabilities is the most 
we can hope for. More importantly for us, we conclude 
from Fig. 2 that it is very difficult to infer the value of 
A G [3.7, 4] from the A-DF of order patterns of length 4. 



B. Order patterns for the skew tent map 



interval P[o,i.2,3] =: (0, 04(A)) is determined by the first 
intersection between f^\x) and /j: (2): 

4 (A) = (6) 

Therefore, the A-DF of tt = [0,1,2,3] (pattern #0) is 
given by 04(A); see Fig. 5(a) for the graphical represen- 
tation of 04(A). The fact that the function 04(A) is bijec- 
tive entails the possibility of estimating A via the relative 
frequency of the order pattern [0, 1, 2, 3]. 

Up to this point it has been assumed that the orbits 
of the various maps considered, were accessible. From a 
more practical point of view, it is also relevant to know 
whether order patterns can be still determined using less 
information about the orbits. This is the case, for in- 
stance, when dealing with the symbolic dynamic asso- 
ciated to a generating partition of the state space. In 
particular, the orbits of maps of T can be transformed 
into binary sequences by the procedure described in [1]. 
In the next section it is explained how to build order 
patterns from those binary sequences. 



The skew tent map, given by 



IV. GRAY CODES AND UNIMODAL MAPS 



fx(x) 



x/X, 

(!-*)/(! 



A), 



if < x < A, 
if A < x < 1, 



(5) 



for x £ [0, 1] and A € (0, 1), belongs to the subclass ^2, 
comprised of those maps of T parameterized by the crit- 
ical point. Furthermore, for the skew tent map fx, the 
maximum value f\(x c ) = fx (A) = 1 is independent from 
A (see Fig. 3). Contrarily to the logistic map, the skew 
tent map does possess a known ergodic invariant measure 
for all A £ (0, 1), namely, the Lebesgue measure on [0, 1]. 
Hence, if P v is given by Eq. (3) with I = [0, 1], the rela- 
tive frequency of the order pattern tt in a typical orbit of 
the skew tent map, coincides with the Lebesgue measure 
of P n , which can be determined analytically. The easiest 
case corresponds to the order pattern n = [0, 1, ... , L— 1] , 
since then P n is an open interval whose left endpoint is 
and whose right endpoint is the leftmost intersection be- 
tween f^ L and /{ L . The relative frequencies of the 
order patterns of length 4, numbered according to Table 
I, are depicted in Fig. 4. In particular, the length of the 



Symbolic dynamics has been thoroughly studied in the 
context of unimodal maps since the seminal contribution 
of Metropolis et al. in [1]. In [3] Gray codes were used as 
a more intuitive way of understanding and applying the 
ideas of [1]. The connection between both approaches can 
be mathematically established with the aid of results in 
[1, 2, 9], as pointed out in [10]. In this section, we address 
the ordinal structure of Gray codes. 

For a unimodal map / defined on the interval / = 
[a, b], any finite orbit {f n (x) : < n < N - 1} 
can be transformed into a binary sequence Gpf(f,x) = 
9(f°(x)) giPix)) . . .g(f N ^ 1 {x)), where g is the step 
function 



if x < x c , 

1 if x > x c . 



(7) 



As x increases from the left endpoint a to the right end- 
point b, the interval / can be partitioned into 2 N subin- 
tervals n , 1 < j < 2 N , each subintcrval containing 
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FIG. 3: The first four iterations of f(x) and the corresponding order patterns of length 4 for the skew tent map, i.e., f^\x) 
for k = 0,1,2,3. 




Order patterns 



FIG. 4: Relative frequencies of the order patterns of length 
L = 4 realized by the skew tent map. 



those x <E I whose orbits have resulted into a given bi- 
nary sequence Gn{I,x). That is, (i) n N ' n /• = for 



f U/f ! 



U £ , and (iii) the binary 



J + k, (ii) / : 

sequences Gjv (/, x) obtained for each x G I 
same. Moreover, the sequences Gjv(/, x\) for xi £ /j 
and G N (f,x 2 ) for x 2 e /j^, 1 < j < 2 N - 1, differ only 
in one bit. Therefore, if we label the 2 N subintcrvals 
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FIG. 5: Order pattern frequency for the skew tent map and 
L = 4 (a) order pattern #0; (b) order pattern #1; (c) order 
pattern #2; (d) order pattern #3. 



Ij N ^ with the 2 W sequences Gjv(/, a;), then the labels of 
continuous subintervals will have only one bit flipped. 

For the sake of illustration, let us consider the skew 
tent map with A = 0.5. In Fig. 6, the division of / = 



[0, 11 into the subintcrvals I, 



(JV) 



each labeled with the 



corresponding binary sequence of length TV, is shown for 
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Rank 


Binary code 


Gray code 





000 


000 


1 


001 


001 


2 


010 


Oil 


3 


Oil 


010 


4 


100 


110 


5 


101 


111 


6 


110 


101 


7 


111 


100 



TABLE II: Correspondence between Gray codes and binary 
codes for three bits. 



TV = 1.2,3. The separation points of the subintervals 

i 



I)" ' are the solutions of the equations 



A/2 2' 



1 < n < N. 



(8) 



If, furthermore, Qn is the set of all binary sequences 
of length ./V produced by a map / g then it is pos- 
sible to endow Qn with a linear order as follows. Given 
Gjv(/, x\) 7^ Gn(J, x 2 ), let i be the first index such that 
g{f l {xi)) 7^ g(f l (x 2 )). Depending on the value of i, we 
distinguish three cases: 

- If i = then Gjy(f,xi) < Gn(J, £2) if and only if 

g{x{) < g(x 2 ). 

- If i > and Gi(f,xi) = Gi(f,x 2 ) contains an even 

number of l's, then Gat(/, x\) < Gn{J,X2) if and 
only if ^(n)) <g(f(x 2 )). 

- If i > and Gj(/, aq) = Gi(f,x 2 ) contains an odd 

number of l's, then Gjv(/, aq) < G]y(f,x 2 ) if and 
onlyifg(f i (x 1 ))>g(f i (x 2 )). 

Gray codes are well known in the context of commu- 
nication theory. The Gray codes of length 3 are shown 
in Table II. The main characteristic of the Gray codes is 
that two consecutive codes differ in only one bit. More- 
over, the order of Gray codes is equivalent to the order 
in Gn (check Table II for N — 3). As a consequence, any 
binary sequence Gjv(/, x) can be interpreted as a Gray 
code of length N [3] , and will be called a Gray code here- 
after. Finally, the order of the Gray codes derived from 
any unimodal map belonging to T is directly linked to 
the order in R of the points x € I. Indeed, it is proven 
in [9, Lemma 4.1] that Gjv(/, x\) < Gjv(/, x 2 ) for some 
N > 1, implies xi < x 2 . This is illustrated in Fig. 6. 



V. GRAY CODES AND ORDER PATTERNS 
FOR UNIMODAL MAPS 

In this section the analysis focuses on the parametric 
unimodal maps of the subclasses T\ or T 2 . In section III 
we elaborated on the dependence of the order patterns 
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FIG. 6: Symbolic intervals for different iterations of the skew 
tent map for A = 0.5. 




FIG. 7: Error rate for the pattern assignment based on Gray 
codes with respect to the one based on the orbit of the logistic 
map. The length of order patterns is L = 4, the length of the 
considered Gray codes is N = 100 and the number of samples 
is 10000. The perfect estimation of the PDF of the order 
patterns of the logistic map is possible for those values of A 
leading to aperiodic binary sequences or to binary sequences 
with period larger than 4, i.e., the length of the considered 
order patterns. 



allowed for those maps with respect to the control pa- 
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rameter. Specifically, we estimated the probabilities of 
order 4-patterns by their relative frequencies in orbits of 
the logistic map (Fig. 2) and of the skewed tent map (see 
Fig. 4) with different parameter settings. Our next goal 
is to reproduce the same dependencies not from the exact 
values of the orbit point ( "sharp orbit" ) , but from the bi- 
nary sequence built as explained in the previous section 
("coarse-grained orbit"). As discussed in that section, 
the definition domain I of / G T splits in 2 N subinter- 
vals when Gray codes of length N are considered. We 
show next that the order patterns of / can also be ob- 
tained comparing Gray codes obtained from its orbits. 

Let G M (f, x) = .9o.9i ■ ■ • 9m-i, 9% € {0, 1}, be the Gray 
code of length M of x G /. Since the Gray codes, 
together with the points x G /, are linearly ordered 
and, moreover, their order relations are equivalent (i.e., 
Goo{f,x{) < Goo(/, X2) iff x\ < X2), we can expect to 
obtain useful information about the order patterns real- 
ized by the sharp orbit Of(x) from the order patterns 
realized by the coarse-grained orbit Gm(J,x), M > 2. 
The procedure is as follows. 

1. Divide the Gray code of length M, Gm(J,x), into 
M - N + 1 Gray codes of length N < M us- 
ing a sliding window of length N. Thus, the 
first Gray code derived from Gm (/, x) is G° = 
.90.91 . . . gN-i = Gjsr(f,x), the second Gray code 
is G 1 = gig 2 -.-9N = G N (f,f(x)), and 
the (M — N + l)-th Gray code is G M ~ N = 
9M—N9M—N+1 ■ ■ ■ .9a/-i = Gn(I, f ~ (x)). 

2. For i = 0, 1, M - N - L + 1, build groups of 
L consecutive Gray codes G l G l+1 ... G i+L ~ l . The 
j-th group defines then the order L-pattern tt = 

7r(£) = [7To,7Tl, . . . ,7Ti_l] if 



fact, remember that the estimation of the probability 
/x(P w ) by the relative frequency of tt G Sl in finite orbits 
of a ^-preserving map, hinges on the ergodic theorem. 
If, furthermore, the convergence of relative frequencies 
to probabilities in the orbits of an ergodic map with re- 
spect to is very slow, a good estimation would require 
exceedingly long sequences — this is what we mean by 
"poor ergodicity" . These errors are shown in Figs. 7 and 
8 for the logistic and the skew tent maps, respectively, 
with 7r = [0,1,2,3], M = 10104, and N = 100. In the 
first case, the value of A lies within the period-3 window 
of the logistic map. In the second case, poor ergodic- 
ity is expected for values of A close to and 1. The 
asymmetry in the error distribution is due to the fact 
that for A ~ 1, the tent map looks like the identity in 
most of / = [0, 1], hence ^0.1,2.3] covers most of /. This 
makes [0,1,2,3] to be the most frequent order 4-pattern 
even when its frequency is calculated using Gray codes. 
Comparison of Figs. 5 and 9 illustrates the accuracy 
of the Gray code-based method for the first four order 
4-patterns (see Table II) of the skew tent map . 
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The order patterns derived using Gray codes need not 
have, in general, similar A-DFs to those derived from 
the sharp orbits. Indeed, order patterns defined by Gray 
codes of length N are built upon the comparison of subin- 
tervals C I (see Sect. IV), rather than comparing 

points of /. The width of the intervals ij^ decreases 
as the length N of the sliding window increases in such 
a way that when N — ■> 00, each one of those intervals 
converges to a single real number. As a result, the error 
in the calculation of the order patterns from Gray codes 
is expected to reduce as N increases. In the context of 
finite-precision computation, the minimum value of N 
necessary to get a reliable approximation of the A-DF of 
an order pattern is related to the precision of the arith- 
metic used. Again, this quantization error decreases as 
TV increases and, consequently, a large value of N may be 
necessary to assure a good approximation of the A-DF. 

Another source of divergences between A-DFs and their 
numerical estimation via finite-length Gray codes maybe 
non-ergodicity or even poor ergodicity. As a matter of 



FIG. 8: Error rate for the pattern assignment based on Gray 
codes with respect to the one based on the orbit of the skew 
tent map. The length of order patterns is L — 4, the length 
of the considered Gray codes is N = 100 and the number 
of samples is 10000. A value of the control parameter above 
0.2 guarantees a perfect estimation of the PDF of the order 
patterns of the skew tent map. 



VI. ESTIMATION OF THE CONTROL 
PARAMETER FOR UNIMODAL MAPS WITH 
CRITICAL POINT DEPENDING ON THE 
CONTROL PARAMETER 

The main characteristic of maps in T2 is that the 
control parameter A determines the value of the critical 
point. Furthermore, from our discussion above, we ex- 
pect that the relation between the control parameter and 
the allowed order patterns of the corresponding dynam- 
ics is specially simple for the pattern n = [0, 1, L — 1]. 
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FIG. 9: Relative frequency of order patterns of the skew tent 
map using Gray codes, when L = 4, N = 100 and the se- 
quences are 10104-bit long: (a) order pattern #0; (b) order 
pattern #1; (c) order pattern #2; (d) order pattern #3. 



Clearly, if the A-DF of this pattern is 1-to-l, then A can 
be pinpointed from that distribution function; otherwise, 
the possible values of A can be reduced to a few candi- 
dates, what can be also acceptable in applications like 
cryptanalysis. In turn, A-DFs can be approximated via 
Gray codes, without previous knowledge of the critical 
point of the map. The bottom line is that the con- 
trol parameter of a map in T2 can be estimated from 
their coarse-grained orbits (in form of Gray codes). The 
specifics depend on the map. 

As an example, consider the skew tent map again. For 
this map, the interval F\o,i,...,l— 1]; i- c -i the set of points 
x € [0, 1] of type [0, 1, ... , L— 1], is determined by the left- 
most intersection of the iterates f^ 2 ail d f\ 1 t where 



KM 



x/X n , if < x < A", 

^ (A™- 1 - .t)/A"- 1 (1 - A), if A" < x < A"" 1 . 

(9) 

Hence P[o,i,...,£-i] = [0, L (A)], with 



vL-2 



2- A 



(10) 



Since this function is 1-to-l in the interval < A < 1 for 
L>2, with ^,(0) = 1/2, L > 3 (O) = 0, and <f> L > 2 (l) = 1, 
it allows to estimate A by estimating 4>l{X) — the length 



of -P[o,i,. ,,,L—i]- Now, from the equation 

A L-3 



2 [2(Z — 2)—{L — 3)A] 



(2 -A) 

[ 0, if A = 0, 
L-l, if A = 1. 



(11) 



it follows that </>£ (A) is a U-convex function on < A < 1 
for L > 2, that converges to on < A < 1 as L — > 00. 
Therefore, the higher L the worse ^l(A) discriminates 
different values of A. Consequently, L — 3, 4 are the best 
choices for a quality estimation of A. 

On the other hand, the ergodicity of the skew tent map 
permits to estimate the length of -P[o,i....,l-i] by estimat- 
ing the relative frequency of the 7r = [0, 1, . . . , L — 1] in a 
typical sharp orbit of the map — or, as we intent, in a typ- 
ical coarse-grained orbit. In the latter case, the choice for 
the parameter N, the width of the sliding window down 
the Gray codes (Sect. V), must be also analyzed. The 
minimum value of N to get a good reconstruction of the 
A-DF of the order patterns, N m i n , depends on the pre- 
cision of the arithmetic used, but it also depends on the 
Lyapunov exponent of the map. If floating point double- 
precision arithmetic is implemented, then N min can be 
determined as function of A by comparing pairs of sym- 
bolic sequences generated from the same initial condition 
and control parameters Ai and A2 such that | A2 — Ai| 
equals the spacing of floating point numbers. As it is 
shown in Fig. 10, the value of N m i n increases with the 
Lyapunov exponent for the skewed tent map. 

Summing up, the estimation of the control parameter 
A € (0, 1) of the skew tent map /a, Eq. (5), can be done 
by counting and normalizing the occurrences of the order 
pattern [0, 1, L — 1], ideally for L = 3 or 4, in a sta- 
tistically significant sample of orbit segments of f\. This 
follows from the following properties: (i) fx is ergodic for 
all A, and (ii) the /^-invariant measure of P[o,i,...,l-i] (in 
this case, the length of the interval P[o,i,...,l-i]) depends 
bijectively on A. In a practical context though, finite 
precision machines are used, and this entails, in general, 
numerical degradation, this meaning that the computed 
orbits, whether of chaotic or non-chaotic maps, depart 
from the real ones. In the case of a very long orbit of 
a chaotic map, the deviation of the numerical simulation 
(locally measured by the Lyapunov exponent of the map) 
will be severe; in such cases, it is preferable to have many 
shorter orbits instead. Even worse, all orbits computed 
with finite precision are eventually periodic. This distor- 
tion of the dynamics, due to finite numerical precision 
and dependence on initial conditions, implies the general 
impossibility of obtaining orbits and invariant measures 
in an exact way. As a matter of fact, all this carries over 
to symbolic dynamics. 

To verify this issue in the case of coarse-grained orbits, 
a sample of Gray codes of the skew tent map, each one 
with the same length but with a different initial condi- 
tion, was generated for every value of A. The underlying 
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(a) 




FIG. 10: Dependence of the width of the sliding window with 
respect to the rate of divergence: (a) Minimum width of the 
sliding window necessary for the reconstruction of the PDF 
of the order patterns from the symbolic sequences of the skew 
tent map; (b) Lyapunov exponent of the skew tent map. 
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FIG. 11: Average deviation and standard deviation in the 
estimation of the PDF of the order pattern #0 for the skew 
tent map with respect to A 2 /(2 — A). 



comprises two stages: 

1. An estimation of A is performed by dividing the 
given Gray code, {ffi}^ 1 , gi e {0,1}, into a 
large enough set of disjoint subsequences of length 
N > 4, say, {gk-N+i}^ 1 for k = 0,1,..., K = 
[AI — N\ — 1. For each such binary subsequence, 
a value of 04(A) is then computed as the relative 
frequency of the order pattern [0,1,2,3] using, of 
course, the Gray ordering (Sect. IV). Let x be 
the mean value of the resulting 04 (A) 's. From 
04 (A) = A 2 / (2 - A), Eq. (6), it follows that the 
control parameter is estimated as 



-x + v x 2 + 8x 



(12) 



sharp orbits were computed with double precision float- 
ing point arithmetic. From this sample of Gray codes, 
the corresponding A-DFs of the order patterns of length 
L = 4 were obtained. The A-DF of the order pattern 
[0, 1, 2, 3] (#0 for short) was calculated as the mean value 
of the A-DFs obtained from the various initial condi- 
tions. This average value is compared to the exact A-DF, 
04(A) = A 2 / (2 — A), in Fig. 11, along with the correspond- 
ing standard deviation. 

Fig. 11 spells out that, in the context of finite precision 
computation, the perfect recovery of the control param- 
eter value using the A-DF of order pattern #0, is not 
feasible in general if one can only resort to Gray codes. 
However, it is possible to locate A up to an uncertainty 
interval. The width of this interval can be upper bounded 
by the standard deviation of the A-DF of the order pat- 
tern #0 since, according to Fig. 11, it is bigger than the 
average error in the estimation of 04(A) for every value 
of A . Therefore, the estimation of the control parameter 



2. If a is the standard deviation of the 04(A) sampling, 
then 



A e (04 (a; - ct), 4 1 (x + a)). 



(13) 



The specifics of this procedure refers to the skew tent 
map, but the general strategy is the same, once a bijective 
A-DF of an order pattern is exactly known. 

In order to establish the accuracy of the procedure, 
some numerical simulations with the skew tent map were 
done. For every value of the control parameter, a group 
of 100 different initial conditions were used in the gener- 
ation of the corresponding Gray codes. For each of these 
binary sequences, the control parameter A was estimated 
as just explained. The mean error of the estimation is 
shown in Fig. 12. The average error lies always above 
10~ 4 , and can only be reduced by the implementation 
of the procedure with extended-precision arithmetic li- 
braries. To prove this claim, the case of the symmetric 
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i 

FIG. 12: Mean error value in the estimation of the control 
parameter of the skew tent map. 

tent map, i.e., the skew tent map for A = 1/2, will be 
now considered. For the symmetric tent map the arith- 
metic is exact. Indeed, if Q.xqX\...xm , %i € {0, 1}, is the 
expansion to base 2 of x <E [0, 1], i.e., 



the shortcomings of finite precision arithmetic and finite 
statistical sampling, but are not inherent to the method 
— as proved with the symmetric tent map. 

10"' f 1 1 , a 




10 L »- 

10 10 10 10 10 

M 

FIG. 13: Dependency of the error in the estimation of the 
rate of occurrences of the order pattern #0 with respect to 
the length of the orbits. 



10 2 



M 

n=0 

(numbers with finite binary expansions are called dyadic 
rationals), then the action of the symmetric tent map 
amounts to a zero-bit dependent left shift, to wit: 

fi/ 2 {0.x xi...x n ...x M ) = 

_ | O.X!X 2 . . . X n +1 ■ ■ ■ X M -lXM, if X = 0, 
I 0.xlx2...x; i+1 ...x* M _ 1 x M , i£x = 1, 

where £* = 1 — x n . Therefore, if x G [0, 1] is represented 
with M bits and xm = 1, the orbit of x collapses to 
after M iterations of /1/2, so M can be considered the 
effective length of the orbits to be used in an estimation 
of A = 1/2. For L = 3, the relative frequency of the or- 
der pattern #0 ([0,1, 2] in this case) was determined for 
a large set of random initial conditions x and increasing 
orbit lengths M. The convergence in average of this rel- 
ative frequency to (fo(l/2) = 1/3 (sec Eq. (10)) as M 
increases, is confirmed by Fig. 13. At the same time, the 
variance of the estimation steadily reduces with M, as 
shown in Fig. 14. In other words, a higher precision of the 
arithmetic used in orbit generation and greater samples 
for the subsequent control parameter estimation, clearly 
improves the results. 

We conclude that the inaccuracies exposed above in 
our method to recover the control parameter of maps of 
Tii based on the order patterns of their coarse-grained 
orbits (specifically in form of Gray codes), are due to 
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FIG. 14: Variance of the error in the estimation of the rate 
of occurrences of the order pattern #0 with respect to the 
length of the orbits. 



VII. CONCLUSIONS 

In this paper it was shown how to rebuild the A-DFs 
of order patterns from Gray codes, the scope being the 
estimation of the parameter A. Gray codes are the 0-1 
sequences that result from the symbolic dynamic of uni- 
modal maps with respect to the left-right partition of the 
state space introduced by the critical point. We have an- 
alyzed the A-DFs of the order patterns of two unimodal 
parametric maps: the logistic map (as representative of 
the subclass T\) and the skew tent map (as represen- 
tative of the subclass T-x). In the case of the logistic 
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map, it turns out that this technique can hardly deliver, 
on account of the complex and many-to-one relation be- 
tween A and those A-DFs. On the contrary, this relation- 
ship is simple, one-to-one, and analytically known for 
ir = [0, 1, ...,L — 1] in the case of the skew tent map. Our 
method improves previous proposals for parameter esti- 
mation of unimodal maps in that a knowledge of the crit- 
ical point value is not needed. However, it demands high 
computational precision and large amounts of data; in 
this regard, we recommend the use of extended-precision 
libraries for good estimations. In the ideal case of arbi- 
trarily high precision, the estimated value of the control 
parameter is arbitrarily close to the real one. 
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